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ABSTRACT 

Elastic properties of the solid regions of neutron star crusts and white dwarfs play an 
important role in theories of stellar oscillations. Matter in compact stars is presumably 
polycrystalline and, since the elastic properties of single crystals of such matter are 
very anisotropic, it is necessary to relate elastic properties of the polycrystal to those 
of a single crystal. We calculate the effective shear modulus of polycrystalline matter 
with randomly oriented crystallites using a self-consistent theory that has been very 
successful in applications to terrestrial materials and show that previous calculations 
overestimate the shear modulus by approximately 28%. 
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1 INTRODUCTION 

Vibrations of neutron stars and white dwarfs are an im¬ 
portant source of information about matter in their interi¬ 
ors. In recent years the elastic properties of neutron stars 
have attracted particular attention because vibrations of 
the crust could provide a mechanism for qu asiperiodic os¬ 
cillations observed in X-r ay burst sources llDunca ni l 19981: 
IStrohmaver fc Wattj[200^ ). The elastic properties of single 
crystals of dense m atter have been calculated in a variet y 
of wo r ks (see, e.g ., [Haensel. Potekhin fc YakovlevI (l2007l) : 
iBaikd (l201ll . [2012 )). However, because of initial variations 
in the local temperature and composition, as well as the 
presence of gravitational and magnetic fields, it appears un¬ 
likely that the solid part of a star is one giant single crystal, 
and the standard assumption in recent work is that matter 
is polycry stalline, with a random d i stribution o f cryst al ori¬ 
entations JOgata fc Ichimar^liggol : lBaikdf201ll . [201^ . The 
question we address here is how to calculate the elastic prop¬ 
erties of such a polycrystalline medium from those of a single 
crystal. 

In white dwarfs and in the crusts of neutron stars at 
zero temperature matter consists of a lattice of nuclei with a 
rather uniform background of electrons, which ensure charge 
neutrality. In the inner crust of neutron stars, the lattice is 
also permeated by a sea of neutrons, but these have little 
influence on the shear properties, which are determined pri¬ 
marily by the Coulomb interaction. 

* E-mail: kobyakov@theorie.ikp.physik.tu-darmstadt.de (DK); 
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In the presence of a uniform background of electrons, 
the lowest energy lattice of nuclei with Coulomb interac¬ 
tions has a body-centred-cubic (bcc) structure and its elas¬ 
tic properties are conveniently described in terms of elastic 
constants cn, ci2, and C44 in the standard Voigt notatiorQ 
or alternatively in terms of the quantities B = (cn -|-2ci2)/3 
(the bulk modulus), C44, and c' = (cn — cyi)I2, which de¬ 
scribe response to three orthogonal distortions of the crys¬ 
tal. The rigidity of the solid to a shear perpendicular to one 
of the crystal axes is described by C44 and that to a shear 
perpendicular to the (1,1,0) and equivalent directions by c'. 
The Coulomb crystal has v ery anisotro pic elastic properties, 
since the anisotropy ratio dZene 3H) 


which is unity for an isotropic medium, is approximately 
7.5 when screening by electrons is neglected. Therefore, it 
is important to know how to relate the elastic properties of 
the polycrystalline medium, which behaves as an isotropic 
medium at long wavelengths and is therefore described by 
just two elastic constants, the bulk modulus and an effective 
shear modulus /ieff- 

The earliest estimate of the effective shear modulus of 


^ For matter at nonzero pressure, which is the case in compact 
stars, it is convenient to define elastic constants in terms of varia¬ 
tions of the Gibbs free energy rather than the Helmholtz free en¬ 
ergy since then the formalism at nonzero pressure becomes iden- 
tic al to the standard t r eatme nt for zero pressure, as pointed out 
bv iMarcus. Ma fc Qiul (l2002h . 
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matter in the crust of a neutron star was ba sed on assuming 
it to be equal to C4d llBavm fc Pinejll97lll . More recently, 
I Ogata fc Ichimanil lll990l l calculated u^fr for polycrystalline 
matter with randomly oriented crystallites on the assump¬ 
tion that it is given by the average of the shear stiffness over 
all possible wavevectors and for polarisation vectors perpen¬ 
dicular to the wave vector, and this prescription has been 
employed in su bsequent w o rk in which many-body effects 
are calculated (lBaikd[201ll . l20f^ f. For a cubic crystal, one 
finds 


2,3 

/ieff.V = -^C + -^C44, 


( 2 ) 


where the subscript V is cho sen si n ce thi s way of averag¬ 
ing wa s proposed earlier by IVojgd (Il887|), a fact pointed 
out by I Johnson-McDaniel fc OwenI (120121 ') in the context of 
the shear properties of mixed phases of quark matter and 
hadronic matter at supernuclear densities. 


2 EFFECTIVE SHEAR MODULUS 

The elastic properties of polycrystalline terrestrial materi¬ 
als is a subject of considerable practical interest and there 
exists an exten sive literatur e on the subject. We now give a 
brief overview. iReussI (Il929l j proposed that, instead of aver¬ 
aging the shear rigidity, one take its harmonic mean. This 
corresponds to taking the shear compliance of the polycrys¬ 
tal, the inverse of the shear rigidity, to be the average over 
wavevectors and polarisations of the shear compliance of the 
single crystal. For a cubic crystal this leads to 

/"S 1 3 1 \-i 

=(-- + -—) . (3) 

For the case of a Coulomb crystal, the Voigt and Reuss pre¬ 
scriptions differ by a factor of about 2.5, so it is im por¬ 
tant t o narrow down the uncertainty in calculating u„ff. iHilll 
(Il952ll showed that, with the assumption that the crystal¬ 
lites have random orientations, the Voigt average gives an 
upper bound to fiefc and the Reuss average a lower bound. 
He proposed that for a more realistic estimate one should 
use either the arithmetic mean of them, which is given by 

neS.H ~ -(/Teff.V + , r ) , ( 4 ) 

which has been commonly used in subsequent wor k, or the 
geome t ric mean iuea .v ■ Follow ing that, iHershevI 
(Il954h . [Kroner! (Il958l ~) and lEshelbvl (Il96ll ~l developed a model 
in which one calculates the strain field around an inclusion 
(representing a single crystallite) in a homogeneous medium 
and then, by averaging over the possible orientations of 
the crystallite, determines the effective shear modulus self- 
consistently (for a pedagogical account of the theory, see 
IdeWitI ll2008h i. In the self-consistent theory one finds 

SfTeff+ (9S + 4c')/reff —3(B-|-4c')c44/refr —6 Rc'c44 = 0.(5) 

The predictions of the self-consistent theory are very robust, 
since very similar values of the effective shear modulus are 
obtained from a rather different approach in which one treats 
the medium as an elastic mediu m with smoo th spatial vari¬ 
ations of the elastic properties ([Krone 3H). 


2.1 Comparison with experiment 


The self-consistent theory gives a rather good account of ex¬ 
perimental measurements of elastic properties of terrestrial 
materials. The anisotropy of matter at high density is ex¬ 
treme, and there are rather few materials with such high 
values. We shall restrict ourselves to cubic materials. 

One system that has been studied extensively is 
the (5-phase of plutonium fl for which A « 7.03 
dhedbetter fc Moment! !l97dj . The experim entally deter¬ 
mine d shear elastic constant, 16.3 GPa ([Migliori et ah! 
[200611 ■ agrees well with the result of the self-consistent 
theory for the single- c rysta l elastic constants measured by 
[Ledbetter fc Moment! ([i976ll . which is 16.2 GPa. Two other 
materials with significant elastic anisotropy that have been 
studi ed in detail are copper, for which A ~ 3.14 ( Ledbett^ 
[l980ll . and stainless steel, with A ~ 3.78 |[Ledbetter|[l984[ l: 
in both these cases, the self-consistent theory gives results 
that agree with experiment to within a few percent. 

The alkali metals also have high elastic anisotropy and 
bcc structure. However, these substances are difficult to 
work with and we have been unable to find recent measure¬ 
ments of the elastic constants of polycrystalline alkali met¬ 
als with well characterised samples. Lithium has A « 8.52, 
and the self -consistent theory predicts fieg = 4.2 GPa 
|[deWit[[2008l i. while the experimental value is 4.0 GPa. The 
behaviour of sodium {A « 6.98) appears to be anomalous: 
the measured shear elastic constant quoted in compilations 
of data is 3.3 GPa, while the self-consistent theory predicts 
2.1 GPa. What is even more surprising is the fact that the 
experimental value is greater than an upper bound in the 
theories, the Voigt average Eq. ((2|, which is 2.8 GPa. With¬ 
out more information about the structure of polycrystalline 
sodium, the reasons for this discrepancy are unclear. 


3 ASTROPHYSICAL DENSE MATTER 


We now apply the results of the previous section to dense 
stellar matter. For illustrative purposes we shall restrict our¬ 
selves to the case of zero temperature and will neglect many- 
body effects. For a Goulomb lattice with a uniform back- 
ground of negative charg e, the elastic constants are given by 


groung or negative cUarg e, 
jFuchj[l93d:fBaiko[[201]M 


c' - 0.09947”^^ ® 

2a 

( 6 ) 

and 


C44 = 0.7424 

2a 

( 7 ) 


^ This phase has a face-centred-cubic structure. However, since 
the theories of average elastic constants depend only on the cubic 
symmetry and not on the detailed arrangement of atoms within 
a cubic cell, this is a relevant system for comparison with theory. 
This phase may be stabilised at room temperature by addition of 
a small amount of gallium, and the results we quote are for this 
dilute alloy. 

^ The elastic co nstants used in the present article are related to 
the quantities in iBaikol ll2Qllf) by c' = (5iiii — S '1122 — -P)/2 and 
C44 = 5i212? where P is the pressure. 
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where a is the length of the side of a cubic cell containing two 
nuclei, njv is the number density of nuclei, Z is the proton 
number of the nucleus, and e is the elementary charge. 

The Voigt and Reuss predictions for the effective shear 


modulus of a polycrystal are therefore 


= 0.4852 

2a 

(8) 

and 



( 9 ) 

while Hill’s suggestion gives 


Meff.g - 0.3462 

2a 

(10) 


We turn now to the self-consistent theory. In contrast to 
the Voigt, Reuss, and Hill expressions, the expression ([5| for 
the shear modulus involves the bulk modulus in addition to 
the two shear elastic constants of the single crystal. In solid 
matter in compact stars, the bulk modulus is dominated 
by that of the elections, UedPe/dne, where rie = Zun is 
the electron density and Pe the electron pressure. This is of 
order for relativistic electrons and is greater than the 

shear elastic moduli by a factor Xj{otZ'^^^') > 10, where 
Q = e^/hc « 1/137 is the fine structure constant. Thus it 
is a good approximation in Eq. to take only the terms 
linear in B. This gives 

3/ieff — C 44 /iefi — 2c C 44 = 0, (11) 

for which the only physically meaningful root (positive ^eff) 
is 

„., = ^(i + ^1 + 24£| (12) 

For the Coulomb lattice, 

/Teff = 0.3778 '"^^^"^" . (13) 

The Voigt average previously used in the astrophysical lit¬ 
erature thus overestimates the shear modulus by approxi¬ 
mately 28%. This effect is much greater than a number of 
many-body effects, such as electron screening. 


4 OBSERVATIONAL CONSIDERATIONS AND 

CONCLUDING REMARKS 

Our results show that, in neutron stars, the method previ¬ 
ously used to calculate the effective shear modulus overes¬ 
timates the frequencies of torsional modes of the crust by 
~ 15% when crystallites have random orientations. 

In white dwarfs, shear waves in the solid will be simi¬ 
larly affected. However, because of the acoustic impedance 
mismatch between the solid and liquid regions of the stars, 
such modes couple weakly to the surface layers of the star, 
and are consequently not easily detectable by obse rvations 
of the stellar surface (iMonteomerv fc Wingedll999ll . 

The self-consistent theory can be applied to quanti- 
ties other than the shear modulus. For the bulk m odulus 
llHershevI 1 19541 and thermal expansion coefficient (iTomd 
Il998h . it predicts for cubic materials that they are the same 
as those of the single crystal. 


The polycrystalline nature of the matter will affect the 
phonon spectrum, since for wavelengths longer than the typ¬ 
ical dimension of a crystallite, the spectrum is that of an 
isotropic solid while for shorter wavelengths it is well approx¬ 
imated by the result for a single crystal, provided the size of 
a crystallite is large compared with the lattice spacing. For 
example, the phonon heat capacity will be proportional to 
T® for temperatures low compared with the Debye tempera¬ 
ture but the coefficient will be different for T > hv±/{LkB) 
and T < Tiv^/[Lks), where L is the typical dimension of a 
crystallite, fcs is the Boltzmann constant, and v± is a typical 
velocity for transverse sound. 

Throughout our discussion we have assumed that crys¬ 
tallites are oriented randomly: further work is required to 
determine whether or not this is a good assumption. This 
requires an investigation of the conditions prior to crystalli¬ 
sation, including temperature inhomogeneities, and nuclear 
composition, as well as the possible role of gravitational and 
magnetic fields. Such studies are also important for deter¬ 
mining the typical dimension of a crystallite. 

A further assumption we have made is that all nuclei 
have the same atomic number Z. If the dispersion in Z is 
large, the assumption that matter may be regarded as a col¬ 
lection of perfect crystallites could be a poor starting point. 

In this Letter we have confined our attention to densi¬ 
ties at which nuclei are spherical to a good approximation. 
In the pasta phases, where nuclei are highly aspherical the 
calculation of elastic properties is more complicated since 
it requires a generalisation of the self-consistent theory to 
liquid-crystal-like systems. 
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